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Abstract: Disturbance by rodents alters the morphologies and nutrients of plants as well as the 
physical-chemical properties of the soils. Changes in plants are considered to be mechanisms of defense 
against the disturbance by rodents. Rodents gnaw on the assimilating branches of Haloxylon ammodendron 
(CA Mey.) Bunge and burrow under the bushes in the desert ecosystems of Xinjiang, China. However, 
eco-physiological responses of different age groups of H. ammodendron to the disturbance by rodents are 
not well understood. In this study, soil physical-chemical properties under the shrubs and the 
above-ground morphological, physiological and biochemical features of assimilating branches of H. 
ammodendron of different age groups (i.e., young, 30-100 cm; middle-aged, 100—200 cm; and mature, >200 
cm) in burrowed and non-burrowed (control) areas were studied in 2018. We found that disturbance by 
rodents significantly increased the crown width and total branching rates of young and middle-aged H. 
ammodendron. Photosynthetic pigment contents of assimilating branches of H. ammodendron were 
significantly reduced under the disturbance by rodents. In term of plant nutrients, the main differences 
among different age groups of H. ammodendron under the disturbance by rodents occurred in the total 
soluble sugar and reducing sugar contents that decreased in young plants, increased in middle-aged plants, 
and did not affect in mature plants. Crude protein and phosphorus contents significantly increased, while 
crude fiber and calcium contents significantly decreased in young plants. Crude fat and calcium contents 
significantly decreased in middle-aged plants. Soil organic matter (SOM), total nitrogen (TN), available 
nitrogen (AN) and available potassium (AK) contents in the topsoil (0—20 cm), which are conducive to 
forming "fertile islands", also increased under the disturbance by rodents. In particular, soil AN and AK 
were the major factors affecting the above-ground morphological characteristics of H. ammodendron in 
burrowed areas. Overall, the response and defense strategies of H. ammodendron to the disturbance by 
rodents differed among different age groups, and the effect of the disturbance by rodents on H. 
ammodendron gradually weakened with the increasing plant age. 
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1 Introduction 


Haloxylon ammodendron (CA Mey.) Bunge, a representative C4 succulent xerophyte grows in the 
desert areas of northwestern China, has xeromorphic characteristics for adaptation to drought and 
high-salinity environments. Recently, researchers found that H. ammodendron provides a high 
nutrient content to herbivorous animals in the desert environments (Kang et al., 2017). The 
assimilating branches of H. ammodendron provide not only nutrients but also water for rodents 
living in the desert environments (Mares et al., 1997). Rodents tend to cut the branches 
horizontally into small pieces and carry them to their burrows for food (Mabry et al., 1977). 
Additionally, fossorial mammals construct burrow systems where food, feces, seeds, plant 
propagules and numerous types of organic material are stored, creating patches rich in organic 
content and affecting the dispersion of soil nutrients (Lacey et al., 2001). Therefore, the life 
activities of rodents are certain to have a profound impact on the growth of H. ammodendron. 

Growth differentiation balance hypothesis projects that plants could progress into a "growth or 
defense" stage, either by investing resources in defenses against herbivores or by out-growing 
competitors (Siemens et al., 2002). Animal feeding damage might change the morphology of 
plants because of the effects of selective feeding on particular plant tissues (Andino and Borghi, 
2018). In certain cases, phenotypic plasticity is regarded as a low-cost defense mechanism that 
has a crucial role in plant—herbivore interactions (Hakes and Cronin, 2011; Kozlov and Zvereva, 
2012; Dong et al., 2017). When the shoot apices of plants have been detached by herbivores, 
damaged plants have been documented to generate axillary branches to improve photosynthesis 
and compensate for tissue damage (Christie et al., 2014; Dong et al., 2017). A study by Zhang et 
al. (2009) also showed that slight damage to H. ammodendron by foraging gerbils promoted an 
increase in the growth of the branches and crown. Nutritional grade of a plant is reflected in its 
probable resistance factor once the nutrient content is outside the optimal array of herbivores 
(Kindler and Robert, 1970; Schweissing and Wilde, 1979). A deficiency in plant nutrients could 
decrease plant palatability (Ngatia et al., 2015), making plants less desirable to herbivores. 
Furthermore, the effects of herbivores on plants might differ among different plant growth stages. 
Studies have shown that disturbance of herbivores increased the defense compounds (e.g., alkaloids) 
of offspring plants (Barton and Koricheva, 2010; Kafle and Wurst, 2019). The increase of defense 
compounds may reduce the investment in growth or the synthesis of nutrients (Ziist et al., 2015; 
Karasov et al., 2017). However, eco-physiological responses of different age groups of H. 
ammodendron to the disturbance by rodents remain unclear. In this paper, morphological 
characteristics of above-ground parts and the level of nutrients in assimilating branches of H. 
ammodendron of different age groups, which is located at the southern edge of the Gurbantunggut 
Desert, under the disturbance of rodents were researched. 

Rodents can modify the physical-chemical properties of soil via their burrowing activities 
(Berke, 2010; Hagenah and Bennett, 2013). Recent studies have shown that the activities of 
burrowing animals play an important role in promoting the formation of "fertile islands", which 
usually manifests as the increases in soil nitrogen, phosphorus and carbon contents (Xu et al., 
2012; Allington and Valone, 2014). Thus, the disturbance by rodents may be beneficial to the 
growth of H. ammodendron. Therefore, we tested the physical-chemical properties of 
below-shrub soil at various depths in burrowed and control areas. Furthermore, the contents of 
leaf photosynthetic pigments are valuable indicators for possible photosynthetic efficiency (Mao 
et al., 2007). Therefore, we also compared the differences in the photosynthetic pigment contents 
of assimilating branches of H. ammodendron of different age groups in burrowed and control 
areas. We hypothesized that: (1) disturbance by rodents would lead to a reduction in 
photosynthetic pigments of assimilating branches of H. ammodendron. To compensate for nutrient 
and photosynthetic losses caused by the disturbance, H. ammodendron might alter its morphology; 
(2) H. ammodendron would weaken rodent feeding on assimilating branches by reducing their 
nutrient contents. However, changes in nutrient contents in assimilating branches of H. 
ammodendron may differ among different age groups; (3) rodent activities would increase the 
contents of organic matter, nitrogen, phosphorus and potassium in the soil under H. ammodendron 
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shrubs, which is likely to stimulate the growth of H. ammodendron to some degree; and (4) 
responses of H. ammodendron to the disturbance by rodents would differ among different age 
groups. In addition, eco-physiological responses of young H. ammodendron to the disturbance by 
rodents would be more dramatic than those of middle-aged and mature H. ammodendron. 


2 Materials and methods 


2.1 Study area 


Study area is located at the southern edge of the Gurbantunggut Desert, at a transitional zone 
between desert and oasis (45°02'N, 85°04’E; 346 m a.s.l.). The study area belongs to a typical 
temperate continental climate with the mean annual precipitation of 117 mm. The mean annual 
evaporation is more than 1980 mm and annual mean temperature is 6.6°C. The maximum 
temperature can reach 43.1°C in summer and the minimum temperature —27.0°C in winter. The 
surface of the landscape is a mosaic of dry riverbed and aeolian sandy landform. In this area, 
xerophytic plants are the most common type of vegetation. The vegetation is dominated by H. 
ammodendron and Haloxylon persicum with approximately 30% coverage (Huang et al., 2015), 
where H. ammodendron is dispersed between the dunes. Rhombomys opimus, Cricetulus 
migratorius, Meriones meridianus and Dipus sagitta are typical rodent species in this area (Xu et 
al., 2012). Amongst them, Rhombomys opimus is the most abundant species. Ma et al. (2018) 
found that the coverage of burrows in the southern edge of the Gurbantunggut Desert in Xinjiang 
was in the range of 18%—20%. 


2.2 Experimental design 


In July 2018, we searched for rodent-burrowed areas along the direction of a longitudinal dune. 
Three 10 mx10 m quadrats were placed in the middle of burrowed area. Taking the burrowed area 
as a center, we set up three 10 mx10 m control quadrats in areas distinct from the edge of 
burrowed area, where there was no (or limited) signs of burrowing, no fresh rodents feces on the 
ground or at the burrowed entrances and no recently plowed sand at the burrowed entrances 
(Tchabovsky et al., 2001). The distances between all quadrats were over 50 m. Two test replicates 
were set at a distance of more than 1 km from the test site, covering three burrowed areas in total 
(control areas). 


2.3 Measurement of morphological characteristics of above-ground plant parts 


Plant height, basal stem diameter, crown width and total branching rate in the quadrats were 
measured for all H. ammodendron. Basal stem diameter was measured approximately 3 cm above 
the ground. Crown width was expressed as the average width of H. ammodendron in the 
east-western and north-southern directions. Total branching rate (Rp) of H. ammodendron was 
calculated using the formula of Whitney (1976) and Steingraeber and Waller (1986). 

Ry = (NT-NS)/(NT-N1), (1) 
where NT is the total number of branches; NS is the highest branch number; and N1 is the first 
branch number. 


2.4 Collection and analysis of assimilating branches 


Five plants from each young, middle-aged and mature H. ammodendron were carefully chosen. A 
total of fifteen plants were selected from burrowed and control areas. H. ammodendron in the 
height range of 30—100 cm was considered as young plant, H. ammodendron of 100-200 cm was 
considered as middle-aged plant, and those over 200 cm was considered as mature plant (Zhang et 
al., 2016; Luo et al., 2017). Assimilating branches of H. ammodendron in all four directions of 
east, south, west and north were collected (annual branches were used to analyze photosynthetic 
pigments, and 1 a to 3 a branches were used to analyze nutrient contents). Assimilating branch 
samples were fully mixed within different areas and age groups, then numbered and brought back 
to the laboratory. We extracted chlorophylls (total, i.e., a and b) and carotenoids with 80% acetone 
and determined their concentrations spectro-photometrically according to the method of 
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Lichtenthaler (1987). Assimilating branch samples were dried in an oven at 65°C until a constant 
weight was achieved, and then crushed by grinder. After grinding, the samples were passed 
through the 40-mesh sieves for further measurement of nutrient content. According to the 
methods of Roe (1955) and Nelson (1944), we quantified total soluble sugar and reduced sugar 
with glucose as the standard. Crude protein content was measured by the Kjeldahl method, and 
the conversion coefficient was 6.25. Crude fat content was determined with diethyl-ether extract 
by the Soxhlet method. Crude fiber content was measured by the weight method (Somogyi, 1952). 
Phosphorous and calcium contents were measured using atomic absorption spectrum-photometry 
(AAS, Type 932GBC, Scientific Equipment Pty Ltd., Australia). 


2.5 Soil collection and analyses 


We carefully selected an area of ground under the canopy cover of fifteen plants each of young, 
middle-aged, and mature H. ammodendron in burrowed and control areas. For each chosen area, 
soil samples were collected at the four compass directions and three depths (0-20, 20-40 and 
40-60 cm) at about 20 cm from the base trunk of the shrub, and samples were mixed to form a 
single composite sample. The soil samples were sealed, labelled and brought back to the 
laboratory for analysis. All soil samples were dried naturally, ground and passed through a 
0.2-mm sieve for analysis. Soil pH value (1:5 soil:water mix) was measured using a pH meter 
(PHS-25 digital explicit, Shanghai Yi Dian Appliance Co., Ltd., China). The following soil 
properties were measured using the methods of Bao (1999): soil organic matter content (SOM, 
KCr207 method), total nitrogen (TN, HClO.z-H2SO, digestion method), total phosphorus (TP, 
Mo-Sb colorimetric method), total potassium (TK, atomic absorption spectrometry), available 
nitrogen (AN, CaCl: extracts analyzed with a flow analyzer), available phosphorus (AP, NaHCO3 
extracts analyzed with the Mo-Sb colorimetric method) and available potassium (AK, NHsOAc 
extracts analyzed by atomic absorption spectrometry). 


2.6 Data analyses 


Softwares of SPSS 19.0 and Origin 9.5 are used for data analysis and plotting. The data in the 
figures are represented by means+SE. The above-ground morphology of H. ammodendron, 
physiological characteristics of assimilating branches, and soil physical-chemical properties in 
burrowed and control areas were compared using an independent-sample T-test method (P<0.05). 
The concentrations of photosynthetic pigments in assimilating branches of H. ammodendron of 
different ages and the physical-chemical properties of the soil at different depths were compared 
by means of a one-way ANOVA (analysis of variance) test and LSD (least significant difference) 
multiple comparisons. Pearson's correlation analysis was used to analyze the correlations among 
the above-ground morphological characteristics. Canoco 5.0 redundancy analysis (RDA) was 
used to examine the associations among plant above-ground morphological characteristics and 
topsoil (0—20 cm) physical-chemical properties in burrowed and control areas (P<0.05). Before 
the analysis, above-ground morphological characteristics and soil physical-chemical properties 
were converted into lg(x+1) form. Detrended correspondence analysis (DCA) was used to analyze 
the above-ground morphological characteristics (Lepš and Smilauer, 2003). The outcomes showed 
that the extended gradients at all axes in burrowed and control areas were 0.325 and 0.328, 
respectively. As these were less than 3.000, they were suitable for linear redundancy analysis. All 
the expansive factors of soil physical-chemical properties were less than 10.000; therefore, they 
were all analyzed by RDA (Braak and Prentice, 1988). Monte Carlo Permutation tests were used 
to progressively choose significant soil physical-chemical properties. 


3 Results 
3.1 Effects of disturbance by rodents on above-ground morphological characteristics of H. 
ammodendron 


Disturbance by rodents reduced the mean height of H. ammodendron, but the difference was not 
significant among different age groups (P<0.05; Fig. 1a). The effect of disturbance by rodents on 
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the above-ground morphological characteristics was more pronounced for young H. 
ammodendron. The base diameter and crown width of young H. ammodendron in burrowed areas 
were significantly larger than those in control areas (P<0.05) by 38.92% and 41.43%, respectively 
(Fig. 1b). Disturbance by rodents significantly increased the crown width of middle-aged H. 
ammodendron by 21.17% (P<0.05; Fig. 1c). Total branching rates of young, middle-aged and 
mature H. ammodendron significantly increased by 20.88%, 38.65% and 48.19%, respectively 
(P<0.05; Fig. 1d) under the disturbance by rodents. 

The height of H. ammodendron showed a close positive correlation with its basal stem diameter; 
nonetheless, this relationship was more significant in the control areas than in the burrowed areas 
(P<0.01; Fig. 1e). A close positive correlation was found between the height and crown width of 
H. ammodendron in control areas (P<0.01; Fig. 1f). However, low-level positive correlation was 
also found between the height and crown width of H. ammodendron in burrowed areas (P<0.05; 
Fig. 1f). It appears that the effect of disturbance by rodents on crown width was greater than that 
on basal stem diameter. 
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Fig. 1 Difference in above-ground morphological characteristic (a-d) of H. ammodendron with different ages 
(heights) in burrowed and control areas. Boxes represent the interquartile range (containing 50% of values), lines 
across boxes represent medians, whiskers indicate the highest and lowest values, and circle indicates outlier. * and 
** indicate significant differences between burrowed and control areas at P<0.05 and P<0.01 levels, respectively. 
(e), Correlation between height and basal stem of H. ammodendron in burrowed (R?) and control (Ro?) areas; and 
(f) correlation between height and crown width of H. ammodendron in burrowed (R?) and control (Ro?) areas. 


3.2 Effects of disturbance by rodents on physiological characteristics of assimilating 
branches 


3.2.1 Nutrient content in assimilating branches 

Effects of disturbance by rodents on the nutrient content of assimilating branches of H. 
ammodendron differed among different age groups (Fig. 2). Disturbance significantly reduced the 
content of most nutrients, i.e., total soluble sugar, reducing sugar, crude fiber and calcium, in 
assimilating branches of young H. ammodendron, and they decreased by 15.21%, 25.78%, 
15.74%, and 17.40%, respectively. However, crude protein and phosphorus contents were 18.04% 
and 15.62% higher in borrowed areas than in control areas (Figs. 2a, c, d, e and f). In contrast, 
total soluble sugar and reducing sugar contents in assimilating branches of middle-aged H. 
ammodendron in burrowed areas were 29.44% and 27.05% greater than those of in control areas, 
whereas, crude fat and calcium contents were 16.36% and 33.01% lower in burrowed areas than 
in control areas (Figs. 2b, c, e and f). For mature H. ammodendron, disturbance by rodents only 
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reduced the phosphorus content in assimilating branches (P<0.05; Fig. 2). Overall, the change in 
nutrient contents in assimilating branches was the lowest in mature H. ammodendron compared 
with young and middle-aged H. ammodendron. 


40 ox 
(a) = Burrowed area Q 
a 32 Œ Control area = 
x y= 
© 4 Crude fiber 2 
5 : S 
= 16 Crude protein = 
Le] es 1 
QO eA v 
8L A gz] 
Z a 
0 LA A d J ] O 
30-100 100-200 >200 30-100 100-200 >200 30-100 100-200 >200 
2.5 
= 2.0 J 
Š z 
5 1.5 £ 
£ 
3 1.0 a 
~ Ò 
9 05 £ 
0.0 
s P 
S4 S 
5 z 
Z3 $ 
E 3 
Fa E 
$! 4 
= o 4 Aa 4 0.0 
30-100 100-200 >200 30-100 100-200 >200 ' 
Height (cm) Height (cm) 


Fig. 2 Difference of nutrient content in assimilating branches of H. ammodendron in burrowed and control areas. 
Bars represent standard error. * and ** indicate significant differences between burrowed and control areas at 
P<0.05 and P<0.01 levels, respectively. 


3.2.2 Photosynthetic pigment content in assimilating branches 

Disturbance by rodents had a significant effect on the photosynthetic pigment content of 
assimilating branches of H. ammodendron. Photosynthetic pigment content of young H. 
ammodendron was significantly reduced (P<0.05; Fig. 3). And in middle-aged H. ammodendron, 
photosynthetic pigment content was significantly decreased, but there was no significant effect on 
carotenoids (P<0.01; Fig. 3a). Chlorophyll a, total chlorophyll (Chlorophyll a and b) and 
carotenoid contents of assimilating branches of mature H. ammodendron significantly decreased 
in burrowed areas (P<0.05; Fig. 3). 
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Fig. 3 Comparison of photosynthetic pigment content of H. ammodendron in burrowed and control areas. 
Different lowercase letters represent significant differences among different age groups (heights) of H. 
ammodendron at P<0.05 level. * and ** indicate significant differences between burrowed and control areas at 
P<0.05 and P<0.01 levels, respectively. 
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3.3 Soil physical-chemical property 


Soil depth had no significant effect on the contents of SOM, TP, AP, TK and AK in burrowed and 
control areas (P>0.05; Fig. 4). Soil TN content was the highest in the surface soil (0—20 cm). 
With increasing soil depth, soil TN content decreased by 55.56% and 25.00% in burrowed and 
control areas, respectively; whereas, soil water content increased by 75.60% and 67.75% in 
burrowed and control areas, respectively (P<0.05; Fig. 4). Disturbance by rodents enhanced the 
"fertile island" effect of H. ammodendron, which mainly manifested as an increase in soil AN 
content at the depth of 0-40 cm, especially in the topsoil (0-20 cm). SOM, TN, AN and AP 
contents significantly increased by 20.09%, 55.57%, 45.06% and 8.94%, respectively, at the 0—20 
cm depth in burrowed areas compared with control areas (P<0.05; Fig. 4). AN and AP contents 
and pH of soil at the depth of 20-40 cm were significantly higher in burrowed areas than in 
control areas (P<0.05; Fig. 4), having increased by 40.39%, 22.71% and 1.57%, respectively. 
However, soil pH at the depth of 40—60 cm in burrowed areas was significantly lower (2.12%) 
than that in control areas (P<0.01; Fig. 4). In addition, disturbance by rodents significantly 
reduced the water content of surface soil (0—20 cm) by 35.56% (P<0.05; Fig. 4). 
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Fig. 4 Difference in soil physical-chemical property of H. ammodendron in burrowed and control areas. 
Different lowercase letters indicate significant difference among different soil depths at P<0.05 level. Boxes 
represent the interquartile range (containing 50% of values), lines across boxes represent medians, whiskers 
indicate the highest and lowest values, and circle indicates outlier. * and ** indicate significant differences 
between burrowed and control areas at P<0.05 and P<0.01 levels, respectively. SOM, soil organic matter; TN, 
total nitrogen; TP, total phosphorus; TK, total potassium; AN, available nitrogen; AP, available phosphorus; AK, 
available potassium. 


3.4 Correlation between soil physical-chemical property and above-ground morphological 
characteristic 


RDA was performed to evaluate the associations between physical-chemical properties of topsoil 
(0-20 cm) and the above-ground morphology of H. ammodendron (Fig. 5; Table 1). In this 
experiment, there was no significant difference in soil water content among the deep soil samples 
(data not shown). Hence, soil water content was not included in the RDA analysis. The first two 
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axes elucidated 70.2% of the morphological characteristics of burrowed areas (Fig. 5a) and 63.8% 
of the morphological features of control areas (Fig. 5b). AN (P=0.002) and AK (P=0.002) had 
important influences on the above-ground morphology of H. ammodendron in burrowed areas and 
the order of influence was AN>AK (Fig. 5a). In control areas, AN (P=0.02), AK (P=0.004), TN 
(P=0.012), and AP (P=0.008) had notable effects on the above-ground morphology of H. 
ammodendron, with the order of influence being AN>AK>TN>AP (Fig. 5b). The results indicate 
that AN and AK in the topsoil are the most factors affecting the above-ground morphology of H. 
ammodendron in the presence or absence of disturbance by rodents. 
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Fig. 5 Correlation between physical-chemical property of topsoil and above-ground morphology of H. 
ammodendron, as determined by redundancy analysis (RDA). (a), burrowed areas; and (b), control areas. A, 
young H. ammodendron; ©, middle-aged H. ammodendron; o, mature H. ammodendron; SOM, soil organic 
matter; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AN, available nitrogen; AP, available 
phosphorus; AK, available potassium. 


Table 1 Forward selection of soil variables during the redundancy analysis 
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Sample plot Soil variable Eigenvalues Explained variance (%) F P 
Burrowed area AN 0.505 50.5 25.5 0.002 
AK 0.309 30.9 11.2 0.002 
TP 0.107 10.7 3.0 0.072 
TK 0.078 7.8 2.1 0.122 
SOM 0.045 4.5 1.2 0.270 
pH 0.038 3.8 1.0 0.354 
AP 0.018 1.8 0.5 0.578 
TN 0.014 1.4 0.3 0.690 
Control area AN 0.426 42.6 18.5 0.002 
AK 0.305 30.5 11.0 0.004 
TN 0.270 27.0 9.2 0.012 
AP 0.256 25.6 8.6 0.008 
SOM 0.101 10.1 2.8 0.104 
TP 0.041 4.1 1.1 0.294 
pH 0.008 0.8 0.2 0.762 
TK 0.002 0.2 <0.1 0.878 


Note: AN, available nitrogen; AK, available potassium; TP, total phosphorus; TK, total potassium; SOM, soil organic matter; AP, 
available phosphorus; TN, total nitrogen. 
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4 Discussion 


4.1 Soil physical-chemical property under the shrub of H. ammodendron 


In this study, rodent activity significantly increased SOM, TN, AN and AP contents in the soil 
(Fig. 4). This is mainly due to the burrow system constructed by rodents under the shrub of H. 
ammodendron, resulting in an accumulation in plant residues, feces and urine in the soil and 
strengthening the "fertile island" effect. However, the difference was most notable in the topsoil 
(0-20 cm) (Fig. 4) because of the range of rodent activity. Studies have shown that gerbils, the 
most abundant rodent species in this area, are mainly active at the 20—30 cm soil depth (Jiang et 
al., 2007; Xu et al., 2012). In addition to their underground activities, the frequent surface 
activities of rodents also lead to an increase in nutrients in the topsoil. The increase in soil AN and 
AP contents at the 20-40 cm depth in burrowed areas may be related to the activity of 
microorganisms. Studies have shown that increase in organic matter and nitrogen contents of the 
soil increases microbial biomass and diversity, improving the bio-availability of minerals (Ewing 
et al., 2007; Jiang et al., 2007). In addition, the activities of rodents are likely to affect the release 
of exudates from the roots of H. ammodendron, thus affecting the physical-chemical properties of 
soil and the biomass and diversity of microorganisms. Moreover, this interaction may be more 
obvious in rhizosphere soil. However, this speculation needs further study. 

As a typical salt-tolerant plant, the salt content of the soil under H. ammodendron is high, due 
to the absorption of water and transport of salt by the shrub's root system ( Munns and Tester, 
2008; He et al., 2018), a phenomenon known as the "salt island" effect (Li et al., 2011). In our 
study, interference by rodents significantly increased the soil pH at the depth of 20—40 cm (Fig. 4). 
The activities of rodents may have an impact on the "salt island" effect under H. ammodendron. 
Damage to the roots of H. ammodendron through the activities of rodents may weaken its ability 
to absorb and transport salt. Our study did not involve the determination of soil salt content, so 
further relevant research is needed. In addition, because the activities of rodents vary among 
different seasons (Liu et al., 2012), the physical-chemical changes in the soil of burrow areas will 
vary depending on the time of soil sample collection. For example, in the season when rodents 
store food, the number of foraging holes will vastly increase, which will affect the 
physical-chemical properties of the soil (Whitford and Kay, 1999). Studies into changes in the soil 
physical-chemical properties due to disturbance by rodents throughout the seasons are needed to 
resolve this question. 


4.2 Above-ground morphological characteristic of H. ammodendron 


We found that AN and AK were positively correlated with the height, base diameter and crown 
width of H. ammodendron and were the main determining factors affecting the growth of H. 
ammodendron in areas disturbed by rodents (Fig. 5a). This might be because the return of 
nutrients to the soil is slow in desert ecosystems (Barthelemy et al., 2018), and N and K sources in 
the rodent feces could be readily absorbed by vascular plants to stimulate plant growth (Sitters et 
al., 2017). However, disturbance by rodents had a minor inhibitory effect on plant height, 
although the difference was not statistically significant (Fig. 1a). The results show that herbivores 
have a small inhibitory effect on the growth of H. ammodendron. Young H. ammodendron in 
burrowed areas had thicker stems than those in control areas (Fig. 1b), and this was probably 
related to the subsequent increase in soil nutrients. However, disturbance by rodents only affected 
the base diameter of young H. ammodendron. In addition, young H. ammodendron disturbed by 
rodents had greater crown widths and total branching rates than those of undisturbed plants (Figs. 
1c and d), which is likely to be related to the increases in the soil nutrient content and 
obtainability. The increased branching rate may also be triggered by amputation, which tends to 
lead to an increase in crown width. Amputation generally leads to the greater branch strength, 
which is induced by axillary meristems released from the apical dominance (Shimizu-Sato et al., 
2009). Furthermore, effect of the disturbance by rodents on the total branching rate of 
middle-aged and mature H. ammodendron was much greater than that of young H. ammodendron 
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(Fig. 1d), which may be due to the fact that rodents favor larger plants when building burrows, 
and they move short distances for small periods of time, staying near burrows most of the time 
(Falkenberg and Clarke, 1998; Rogovin et al., 2003). Hence, larger H. ammodendron display a 
greater branching response to damage. The increased crown width and branching help to increase 
the leaf area of H. ammodendron, allowing it to absorb more sunlight, CO2 and O2. To a degree, 
this is also because of the morphological flexibility of H. ammodendron in response to the 
disturbance by rodents. 


4.3 Photosynthetic pigment and nutrient contents of assimilating branches of H. 
ammodendron 


Disturbance by rodents significantly decreased the photosynthetic pigment content of annual 
assimilating branches of H. ammodendron, except for the carotenoid content in assimilating 
braches of middle-aged H. ammodendron (Fig. 3). It is well known that jasmonic acid (JA) or 
methyl jasmonate (JA-Me) promotes chlorophyll degradation in leaves (Jung, 2004; Li et al., 
2016). Both mechanical injury and herbivore attack may lead to a rapid accumulation of JA or 
JA-Me in plants (Campos-Vargas and Saltveit, 2010; Nabity et al., 2013). Hence, we speculate the 
reduction in chlorophyll content triggered by disturbance by rodents might be associated with 
JA-induced chlorophyll degradation after assimilating branches were fed. Furthermore, 
mechanical damage to chloroplasts and reactive oxygen species-mediated lipid peroxidation of 
chlorophyll pigments may also be the influencing factors of chlorophyll content reduction 
(Sivritepe et al., 2009). However, this phenomenon requires further study. 

Compared with control areas, there were no significant correlations among the content of 
photosynthetic pigments and the total soluble sugar or reducing sugar of assimilating branches in 
burrowed areas (Fig. S1). This indicates that the decrease in photosynthetic pigment content is not 
the main reason for the changes in soluble sugar and reducing sugar contents. In our study, 
disturbance by rodents significantly increased the contents of total soluble sugar and reducing 
sugar, while significantly decreased the content of crude fat in assimilating branches of 
middle-aged H. ammodendron (Figs. 2b, e and f). It may be related to the mechanical damage of 
photosynthetic tissues by herbivores, changing the status and distribution of carbohydrate in 
plants (Koch, 1996). In burrowed areas, total soluble sugar and reducing sugar in assimilating 
branches may be mostly used as osmoregulation substances for middle-aged H. ammodendron. In 
contrast, total soluble sugar and reducing sugar did not increase, but instead significantly reduced, 
in assimilating branches of young H. ammodendron after being disturbed by rodents (Figs. 2e and 
f), possibly because young trees lack the well-developed root and stem food stores that are 
omnipresent in mature plants (Barton, 2016). Hence, after being browsed, young H. 
ammodendron reduced the nutrient content, such as total soluble sugar, reducing sugar, crude fiber 
and calciumin assimilating branches (Fig. 2) to discourage feeding by rodents. According to the 
optimal allocation theory, extremely defensive plants produce chemical substances to prevent 
animal feeding, which might decrease the availability of nutrients for growth and development 
(Züst et al., 2015; Karasov et al., 2017). However, our study only focused on the levels of 
nutrients, and any changes in the levels of chemical defense substances needs further study. 

Modification of nutrient contents in assimilating branches of H. ammodendron due to the 
disturbance by rodents was only significant in young H. ammodendron. Ulappa et al. (2014) 
reported that browsed plants with a long-term occupancy were more likely to have the greater 
crude protein content. In our study, there was a significant positive correlation between crude 
protein and soil N (Table S1). However, disturbance by rodents only significantly increased the 
crude protein content in assimilating branches of young H. ammodendron (Fig. 2a). The nutrient 
contents in assimilating branches of mature H. ammodendron were virtually unaffected, except 
for phosphorus contents (Fig. 2). These findings are comparable to the results reported by 
O'Reilly-Wapstra et al. (2005) and Sotelo et al. (2014). The age at which plants are fed on also 
affects the optimum defense response (Karasov et al., 2017). Mature plants of the same species 
might have a stronger capability to defend themselves than young specimens. Thus, the degree of 
disturbance by rodents seen in burrowed areas of this study was inadequate to have a significant 
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impact on the nutrient composition of assimilating branches of mature H. ammodendron. 

It should be noted that the amount of feeding on H. ammodendron by rodents varies with 
season and burrow density (Liu et al., 2012; Xu et al., 2012). H. ammodendron was much more 
abundant than any other shrub species and had become the main food source for rodents in the 
study area during autumn, when most herbaceous plants withered (Liu et al., 2012). We conducted 
the experiment in summer when the burrow coverage was 18%-—20% in this area (Ma et al., 2018). 
However, with a change in the degree of disturbance, eco-physiological changes in H. 
ammodendron may fluctuate. The results of our study preliminarily examined some of the 
eco-physiological responses of H. ammodendron to the disturbance by rodents in the southern 
edge of the Gurbantunggut Desert. 


5 Conclusions 


Rodent feeding decreased the photosynthetic pigment contents of H. ammodendron but increased 
the volume of branches and crown to cope with the losses of plant tissue, nutrient and 
photosynthetic machinery. In response to the losses caused by rodent activity, young H. 
ammodendron tended to reduce the nutrient content of assimilating branches to discourage 
herbivores from feeding on the branches. The effect of increased soil nutrient on young H. 
ammodendron was more dramatic than those of on middle-aged and mature H. ammodendron. 
Middle-aged H. ammodendron tended to invest more energy in defense, increasing their chemical 
defense substances, such as total soluble sugar and reducing sugar and decreasing the nutrient 
content of assimilating branches. The response of mature H. ammodendron to rodents was the 
least pronounced of the age groups. Overall, H. ammodendron with different age groups has 
varying degrees of tolerance to the disturbance by rodents. Defense system of H. ammodendron 
increased with the growth of plants and the effect of disturbance by rodents on H. ammodendron 
gradually decreased with age. 
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Fig. S1 Pearson's correlation coefficient among physiological indices of assimilating branches of H. 
ammodendron in burrowed and control areas. * and ** indicate significant differences between these indices at 
P<0.05 and P<0.01 levels, respectively. TSS, total soluble sugar; RS, reducing sugar; CF, crude fiber; EE, crude 
fat; CP, crude protein; P, phosphorous; Ca, calcium; Cha, Chlorophyll a; Chb, Chlorophyll b; Car, carotenoid; 
Total Chl, total chlorophyll. 


Table S1 Correlation of nutrient between topsoil and assimilating branches of H. ammodendron in burrowed 
area 


SOM TN TP TK AN AP AK 

TSS 0.58 —0.32 -0.55 —0.21 -0.35 0.50 —0.46 

RS 0.55 —0.38 —0.52 —0.17 —0.40 0.35 —0.29 

CF 0.06 —0.47 0.03 0.26 —0.35 —0.06 0.34 

EE —0.38 -0.82™" 0.12 —0.30 -0.79" —0.43 0.31 
CP 0.11 0.62" 0.70" 0.15 0.75" 0.18 0.78" 

P 0.03 0.61 0.54 0.11 0.78" 0.71" 0.40 
Ca 0.14 -0.74" 0.54 0.43 -0.88"" 0.31 -0.82™" 


Note: TSS, total soluble sugar; RS, reducing sugar; CF, crude fiber; EE, crude fat; CP, crude protein; P, phosphorous; Ca, calcium; SOM, 


soil organic matter; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AN, available nitrogen; AP, available phosphorus; AK, 
available potassium. * and ** indicate significant correlations at P<0.05 and P<0.01 levels, respectively. 


